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F. Erdogan and B. Aksel 


ABSTRACT: In this paper after giving a general description of the 

line spring model it is extended to cover the problem of interaction of mul- 
tiple internal and surface cracks in plates and shells. The shape functions 
for various related crack geometries obtained from the plane strain solution 
and the results of some multiple crack problems are presented. The problems 
considered include coplanar surface cracks on the same or opposite sides of 
a plate, nonsymmetrically located coplanar internal elliotic cracks, and 
in a very limited way the surface and corner cracks in a plate of finite 
width and a surface crack in a cylindrical shell with fixed end. 

KEY WORDS: Stress intensity factor, line spring model, surface crack, inter- 
nal crack, three dimensional crack problem, part-through crack, cracks in 
plates, cracks. in shells. 


1. Introduction 

The analysis of a part-through crack in a component which may locally 
be represented by a "plate" or a "shell" is certainly one of the important 
problems in fracture mechanics. The general problem is one of a three- 
dimensional crack in a solid with bounded geometry where there is a strong 
interaction between the stress field disturbed by the crack and the bounding 
surfaces of the medium. Even under the assumption of linear elasticity a 
neat analytical treatment of the problem seems to be intractable. The 
existing solutions, therefore, rely very heavily on the techniques of com- 
putational mechanics. In this respect the standard technique has been that 
of three-dimensional finite element (see, for example, [1 ]- [5] for some of 
the typical contributions). Other numerical techniques used have been the 
boundary integral equation method [6], and the alternating method (see, for 
example, the article by Shah and Kobayashi in [7]). The finite element and 
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the alternating methods have also been combined in a new hybrid technique to 
treat three dimensional elliptic crack problems [8], [9], 

If one is dealing with a relatively thin-walled structure, containing a 
part-through crack in a plane perpendicular to the bounding surfaces of the 
medium, the three-dimensional crack problem can be made analytically tractable 
under two important approximating assumptions. The first is the representa- 
tion of the structure by a "plate" or a "shell" and the second is the treat- 
ment of the part-through crack by a "line spring model". Through the use of 
a plate or a shell theory the coordinate in the thickness direction is sup- 
pressed and the basically three-dimensional elasticity problem is rendered 
two-dimensional. Modeling of the crack by a line spring, on the other hand, 
not only lends itself to plate or shell treatment but also preserves the 
basic plane strain character of the stress field along the crack front (every- 
where except near the ends). In a relatively thin-walled structure contain- 
ing a part-through crack the net ligament in the plane of the crack would 
generally have a constraining effect on the crack opening displacements. The 
physical concept underlying the line spring model which was first proposed 
in [10] consists of approximating the three dimensional crack problem by a 
coupled membrane-bending problem through reducing the net ligament stresses 
to the neutral axis of the plate or the shell as an (unknown) membrane load 
N and bending moment M. In the resulting problem the crack surface displace- 
ments are also represented by two lumped quantities, namely the crack opening 
displacement <$ and the crack surface rotation e measured at the neutral sur- 
face. The unknowns N, M, 6 and e are now functions of a single variable x-j , 
the coordinate along the crack in the neutral surface (Fig. 1). The comple- 
mentary pairs of functions (N,M) and (<$,e) defined along the crack are not 
independent and are assumed to be related through the corresponding plane 
strain problem for the cracked strip. The functions (N,M) or (S,e) are deter- 
mined from the related mixed boundary value problem for the plate or the 
shell with a through crack in which N and M are treated as unknown crack sur- 
face loads. After determining N and M the stress intensity factor at a 
given location x-j along the crack is calculated from the corresponding two 
dimensional elasticity solution of the cracked strip lying in a plane per- 
pendicular to x-| (Figs, lc, 2, and 3), 
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In the original model classical plate bending theory was used and the 
problem was formulated with N and M as the unknown functions ([10], see also 
Rice's article in [7]). In the application of the model the main problem is 
the solution of a plate or a shell containing a through crack and subjected 
to membrane or bending loads. On the other hand, to have an asymptotic solu- 
tion around the crack tips in a plate under bending which is compatible with 
the elasticity results the necessity of using a higher order plate bending 
theory such as that of Reissner's [11], [12] has now been well-established 
[13]-[17]. Unlike the classical theory (in which one can use one less boun- 
dary condition than needed through the Kirchhoff assumption), the transverse 
shear theories of plates and shells can accommodate all stress and moment 
resultants (or their displacement complements) on the crack surfaces separ- 
ately, that is, one can specify three boundary conditions in plates and five 
in shells. The consequence of this is that the transverse shear theories 
give asymptotic results which are identical to those obtained from plane 
strain and anti plane shear elasticity solutions of the crack problems [16], 
[17], whereas not only the angular distribution of the stresses given by the 
classical plate and shell theories are different than the elasticity results, 
but for skew-symmetric problems even the powers of singularity (in the trans- 
verse shear stress) are not in agreement [16]. One may note that compatibility 
of the asymptotic solutions obtained from the transverse shear theories with 
the elasticity solutions is not restricted to crack problems. As shown in 
a recent study [18] the agreement is valid for general wedge-shaped plates 
of an arbitrary angle. 

An additional change introduced into the line spring model as presented 
in this paper is the use of the "displacements" ( 5 , 0 ) rather than the stress 
and moment resultants (N,M) as the unknown functions in the integral equa- 
tions. This is simply a matter of convenience, as 6 and o turn out to be 
more natural unknowns in formulating the problem. 

After a period of neglect recently there seems to be a renewed interest 
in the applications of the model. The reasons for this are that judiciously 
used the model can give results with acceptable accuracy, requires (compara- 
tively speaking) only nominal computational effort, is highly flexible with 
regard to crack profile, and, as will be indicated in this paper, can be 
routinely extended to treat multiple crack problems of varying geometires . 
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Some recent applications of the model to a variety of plate and shell prob- 
lems may be found in [19]-[27], 

In this paper the line spring model is generalized to cover multiple 
co-planar cracks of arbitrary orientations including internal as well as the 
surface cracks and the related stress intensity shape functions are given. 
The model is applied to plate and shell problems by using Reissner's trans- 
verse shear theory and the results of some typical examples are discussed. 


2. Description of the Line Spring Model 

Let us first consider the problem of a surface crack in a plate or a 
shell under Mode I loading condition (Fig. 1). Referred to the local coordi- 
nate system it will be assumed that u-j , u 2 , u^ are the displacement compo- 
nents, p-j , $2 are the rotations of the normal to the neutral surface (in, 
respectively, x-jX^ and X2X3 planes), and N-j, M^., and , (i,j=l,2) are, 
respectively, the membrane, bending, and transverse shear resultants. Let 

CO OO 

N 2 2 and M22 be external loads applied to the structure away from the crack 
region (Fig. 1) and define 

= N 22 /h , m ro = 6M 22 /h 2 . (1) 

Similarly, let N and M represent membrane and bending resultants statically 
equivalent to the net ligament stress ^(xpOjXg) (Fig. 1c) and define 

a(x 1 ) = N(x-|,0)/h , m(x-|) = 6M(x-|,0)/h 2 . (2) 

Defining now the unknown functions 

f l ^ X 1 ^ = ¥x^j" B 2< x 1 ,+ °) » f 2^ x l^ = sx^ u 2^ x l» + °) » 


and referring to, for example, [24] and [17] for details, the 
problem for the structure under the applied loads (1) and (2) 



k i j(x 1 ,t) f j(t)dt = g i (x 1 ) 

a 

r 


(-a<x-|<a, i=l ,2) , 


f j (t)dt = 0, (0=1 ,2) , 


-a 


through crack 
may be expressed 

(4) 

(5) 


- 4 - 



where 

g-j (x-j ) = [-moo+m(x 1 )]/6E , g^x-j) = [-a 00 +a(x 1 )]/E . (6) 

In shells the integral equations (4) are coupled (i.e., and k^ are non- 
zero), whereas in plates bending-membrane coupling would be through the in- 
put functions m and a only. Note that m and a tend to close the crack sur- 
faces while the external loads m*, and a m tend to open them. 

The first approximating assumption made in developing the line spring 
model is that the crack is a through crack and the constraint caused by the 
net ligament stress a 22 (x.| .O.x^) tending to prevent the crack from opening 
and rotating may be accounted for by applying the membrane and bending 
resultants N and M on the crack surfaces. The second major assumption states 
that the stress intensity factor along the crack front at a location x-| may 
be approximated by the corresponding plane strain value obtained from a strip 
which contains a part-through crack of length L(x-j) and is subjected to uni- 
form tension N(x-|) and bending M(x-j) away from the crack region (Figs. 1 and 
2). The main problem in the development of the model is expressing the 
functions N and M in terms of the unknown functions f^ and f 2 and it is the 
second assumption which makes this possible. 

To obtain N and M in terms of f-| and f 2 we now express the energy avail- 
able for fracture in two alternate forms. First we note that at a location 
x-| along the crack front by using the crack closure concept the energy avail- 
able for fracture may be expressed as 

G = (U-V) = K ] 2 (7) 

where L ( x-j ) is the crack size, is the stress intensity factor, U is the 
work done by the external loads, V is the strain energy, and E and v are the 
elastic constants. From the solution of the plane elasticity problem the 
stress intensity factor for an edge crack shown in Fig. 2 is obtained as 
follows: 

K-, = /FT [og t (s)+mg b (s)] , s = L/h (8) 

where a and m are given by (2). For analytical convenience the shape func- 
tions g t and g b may be expressed as 
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s= L/h . 


(9) 


6 p4 6 *j 

g t (s) = /irs z b.s , g b (s) = /S? z c.s , 

o 0 

(*) 

The calculated values' ' of K-j and the coefficients and c- obtained from 

a least square curve fit are given in Table 1. 

Referring to Fig. 2, we now let d6 and de be the changes in the "load 
line displacements" 6 and e (corresponding to "loads" N and M) as the crack 
length goes from L to L+dL under "fixed load" conditions. From Fig. 2b it 
then follows that 


dU = Nd6 + Mdo (10) 

dV = f[N(<5+ds)+M(e+de)] -f (N6+Me) = f(Nds+Mde) , (11) 

giving the energy available for crack growth dL as 

d(U-V) = f (NdS+Mde) . (12) 


Observing that for constant N and M and varying L 

d6 = ff dL , de = ff dL , 

from (12), (13) and (7) we obtain 

1 / M 3 5 ■ M 90 


Ji-tU-V) . e »^(Nf ♦ Mf) 


(13) 


(14) 


If we now define the matrices 


= Up = = U,) ■ rp . s(s) = (g,p = 


_ r he/6- 


9b g h9t 

Ct 't‘ ] ’ 05) 


substituting from (2) and (8) into (14) we obtain 


G = ( 1 ■ - h)T T Gx = (f) r 


T 9a) 

9L » 


(16) 


giving 

(*) The values of K, shown in Table 1, are from [28] and are considered to 
be accurate. 
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( 17 ) 


3 to _ 2(1 - v 2 ) 

3L " E bT * 


Note that 6 is a function and t is independent of the variable L and w=0 for 
L=0. Thus, from (17) it is seen that 

L 


u> 


- 2(W) , 

E V 


GdL) - 


(18) 


f^ (t)dt 


(19) 


Referring to (3), (5) and Fig. 2, if we observe that 

X 1 X 1 
6 = 2u 2 (x 1 ,0) = 2 | f 2 (t)dt, 6 = 2B 2 (x-j,0) = 2 
-a -a 

equation (18) gives the desired relationship between the pairs of complemen 
tary quantities (m,a) and (f i *f 2 ) which may be expressed as 

L 


= r m l = ( 

L o J VY? { 


GdL) 


-l 


fe/ 1 f-j (t)dt 

D -a 1 


Z 1 f 2 (t)dt 
— 6 


( 20 ) 


It should also be observed that since the crack depth L(x^) is a known func- 
tion of x.| , the coefficient matrix C defined by 

L 


C = ^ c ij (x l^ = TYF ( 


GdL) 


-l 


(21) 


would also consist of known functions Cij(x-[) of x-j . 

Substituting now from (20) and (21) into (4) and rearranging, we obtain 
a *1 

| [k 11 (x 1 ,t)f 1 (t)+k 12 (x 1 ,t)f 2 (t)]dt c 11 (x 1 )| f-j (t)dt 


-a 


1 , . m 

1T C 12 X 1 f 2 (t)dt = - , (~a<x-|<a) , 

-a 

[k 2 -| (x-j ,t)f-j (t)+k 22 (x-j ,t)f 2 (t)]dt - g- c^ (x-j ) 
x i 

f °CO 

— c 22 (x-| ) I f 2 (t)dt — - -g- > ( _ a<x-|<a) • 


f^ (t)dt 


(22a, b) 


-a 


- 7 - 



After solving (22) for f^ and f^ the stress intensity factor K-|(x) may be 
obtained from (20) and (8). 

If the plate or the shell contains collinear surface cracks in x-|Xg 
plane along a k <x-j<b k , (k=l,...,n), the integral equations (8) remain essen- 
tially unchanged and the problem can be solved by defining the unknown func- 
tions f-j and f 2 given by (3) for each crack separately. 

3. Internal and Multiple Cracks 

The line spring model described in the previous section can be applied 
to any coplanar part-through crack problem provided on any cross-section 
parallel to the XgXg plane there is only one net ligament and one crack, if 
any (Figs. 1c and 2). The model can also be extended to apply to part- 
through cracks involving more than one net ligament as in internal cracks 
(Fig. 3) or more than one crack as in coplanar surface cracks on both sides 
of the plate or the shell. Here the major difficulty lies in the fact that 
in such cases usually there is more than one stress intensity factor that 
must be expressed in terms of more than one dimensionless variable. Thus, 
aside from the generalization of the basic concept, the problem reduces to 
sufficiently accurate parametrization of the stress intensity factors. 

Consider, for example, the nonsymmetric crack geometry shown in Fig. 

3. The integral equations for the corresponding through crack problem is 
again given by (8). The difference between the two problems is in express- 
ing the resultants a(x^) and m(x-j) in terms of the unknown functions f-| and 
f^ defined by (3). At each cross-section we note that there are two dimen- 
sionless length parameters, L(x^)/h and d(x-|)/h. The problem can be, how- 
ever, simplified quite considerably if we restrict the discussion to cracks 
which are symmetric with respect to x 3 =d plane (Fig. 3). Thus, defining 
the plane internal crack by 

L(x-j) L(x.|) 

-a<x-|<a, X2=0> d ^ — < x.j < d + — ^ — > (23) 

in the discussion that follows it will be assumed that d is independent of 
x-|. Referring to Fig. 3b, if and Kg are the stress intensity factors 
at the crack tips A and B obtained from the plane elasticity solution, for 
an increase dL in the crack length the energy increment available for 
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fracture (as determined from crack closure) may be expressed as 


d(U-V) = ^ [K a 2 ~ + Kg 2 f-] , (24) 

or 

G = ^ (U-V) = ^ 0< a 2 +K b 2 ) . (25) 

We now note that as long as there is only one variable L representing 
the crack size, the argument leading to the expression of G in terms of 
(N,M) and ( 6 , 0 ) will remain unchanged (see eqs. (10)-(14) ) and from (14) 
and (25) it follows that 

I < N ft + M If) - IT ' K A 2 +K B Z > • (26) 

The solution of the plane elasticity problem shown in Fig, 3b is available 
[29] and « A and Kg can be expressed in terms of certain shape functions as 
follows: 

K a = 4T [eg t (s) + mg Ab (s)] , 

Kg = [eg Bt (s) + mg Bb (s)] , (27) 


s = L/h, a = N/h, m = 6M/h 2 . 

The shape functions are, in turn, expressed by 

g At (s) = ^ b flj s 21 , g Ab (s) = /S’ 2 c Aj s f , 

n 2 . n . 

g Bt (s) = 4? e b Bi s , g Bb (s) = /Its e c Bi s , s = L/h , 


where for a given d the coefficients are obtained from [29] by using a least 
square curve fit [30], 

For d=0, K A =Kg in tension and K A =-Kg in bending and the corresponding 
membrane and bending stress intensity factors as well as the coefficients 
b. and c- are given in Table 2 [30], In the case of nonsymmetrically located 
internal crack, for six different values of d/h the coefficients of the 
shape functions as defined by (28) are shown in Table 3 [30], 
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By examining the equations (14) through (22), from the derivation given 
in the previous section it may be seen thatthe only difference between the 
models representing the surface crack and the internal crack will be in the 
matrix G(s) defined by (15) for the edge crack. In particular for the inter- 
nal crack problem (20) will remain valid provided the matrix G is evaluated 
from 

, g Ab +g Bb g Ab g At +g Bb g Bt 

G(s) = j 2 (29) 

g Ab g At +g Bb g bt g At +g Bt 

which follows from (26)-(28) and (16). 

Another special case is that of two coplanar surface cracks symmetrically 
located on the opposite sides of the plate (Fig. 4). Going through the 
argument step by step one may easily show that, except for the shape func- 
tions, this case is identical to the internal crack problem with d=0 (Fig. 

3). If the plate is under membrane loading only, because of symmetry no 
bending would take place and there is no need for the bending components of 
the shape functions, g Ab , g Bb . For the corresponding symmetric edge cracks 
of depths L/2 the stress intensity factors and the coefficients b. for the 
membrane shape function are given in Table 4 [30]. 

For more general crack geometries the problem can be rather complicated. 
Consider, for example, the nonsymmetric case of the surface crack problem 
shown in Fig. 4. Let the two cracks be defined by reasonably smooth arbi- 
trary functions L-j(x-j) and ^(x-j) and again designate the crack tips at an 
x-| = constant plane by A and B (corresponding to cracks L-j and L 2 , respec- 
tively). The energy available for incremental crack growths dL^ and dL 2 
may again be expressed in the following alternate forms: 

d(U-V) = ^ [K A 2 dL 1 +K B 2 dL 2 ] , (30) 

d(U-V) = ^[(N ^j- + M ^-)dL 1 + (N -^- + M ^-)dL 2 ] . (31) 

If we now define the matrices as in (15) replacing G by G^ and Gg (where 
the shape functions as defined in (27) would be functions of the variables 
s -j =L-j /h and s 2 =l 2 /h), from (30) and (31) it can be shown that 
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(32) 


3ui 

3L, 


dL l + ^ L dL 2 = da) = (W'V 1 ^' 


Again by observing that io=0 for L-|=0=L2 and L-j =Li ( ) , L 2 =L 2 (xi), from 
(32) we find 


CO 


(X,) , 20^v£i ( 


G A dL l + 


G B dL 2 ) T 


(33) 


From (33) and (18)- (21 ) it then follows that the integral equations (22) are 
still valid provided the matrix C defined by (21) is replaced by 

1-1 1-2 1 

c = ( c ij(X])) = 2(1 h?) G A dL l + G B dL 2^ * ^ 


Of course the main difficulty in problems such as the one described 
above is that they require a complete two-way parametrization of the stress 
intensity factors or the determination of the shape functions g _ (a=(A,B), 

Ot p 

3=(t,b)) as functions of two variables s^L-j/h and s 2 =L 2 /h. 


4. Applications and Some Results 


In order to apply the line spring model to coplanar multiple cracks 
the integral equations need to be cast in a more convenient form. First we 
note that for the case of through cracks in plates and shells derivation of 
the integral equations for colli near multiple cracks is no more difficult 
than for a single crack. In fact, if no symmetry with respect to the coordi- 
nate x^ is required (Fig. 1), the expressions of the kernels for the two 
cases are identical. If, for example, the structure contains P through 
cracks in x-jX^ plane along ap < x-| < hp, p=l,2,...,P, the integral equations 
for a plate or a shallow shell (replacing (4) and (5)) may be expressed as 


z 

p=l 


P 2 


£ k i j ( X 1 ,t)f j(t)dt - 9 ^ ( x i ) * X-j E E (Sp,bp) * 


fj(t)dt = 0, j = 1,2, p = 1,...,P , 


(35) 


(36) 
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where the functions again represent the crack surface loads and in the 
case of part-through cracks the definitions given by (3), (6), (1) and (2) 
are still valid. We note that there are really 2P unknown functions in the 
problem, namely f-j and f2 on each one of the P cracks. Thus, for the pur- 
pose of solving the integral equations, (35) can be written in a more con- 
venient form by defining 

W = f ip( x i ) > a p <x l <b p » P=1 » — »P» 1=1.2, (37) 


m(x-j) = m p (x 1 ), a p< x i <bp, p=l,...,P, 


a( Xl ) = o p (x-,), a p<x-j <bp, p=l,...,P . 


(38) 

(39) 


We also observe that for each crack the net ligament resultants m and a are 
related to f-j and f 2 through equations such as (20) and (21) where the 
matrices G and C are dependent on the local (two dimensional) crack geometry 
and C is a function of x-j. This means that, using (21), (20) may be replaced 
by 


i- x 


m_ _ 

p - c £ ’ EC u,> 


h 

F 


X rl 


1 


f lp (t)dt 


f 2 p (t)dt 


, Sp^x-j <b p * P _ 1.»»*»P« 


(40) 


If we now express the integral equations for each interval (ap.b ) separately, 
from ( 35 )- (40) and (6) we obtain 
b 


P fP 2 

Z 


P=1 


.f k l j^ x l * t ) f jp^ t ) dt ” 36 C 11^ X 1^J f lr^^ dt 
3 — * a 


1 


- g c 12^ x l ^ 


X rl 


m 


f 2r (t)dt = - ef . ( a r< x l < b r * r= l»-**>P) 


(41) 
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( 42 ) 


(43) 

y* 

where the functions c ij(x-|) are determined from (21) or (34) by using the 
geometry and the shape functions for the rth crack. 

It is seen that once the kernels k. . corresponding to the through cracks 

* J 

are determined and the part through crack profiles L r (x-|) (or L r -| and L^) 
are specified, the integral equations (41) and (42) (subject to conditions 
(43)) may be solved for f. (t) and equations (40) and (8) or (27) would 

J ' 

then give the stress intensity factors as functions of x-,. The kernels k- . 

• ' sJ 

for various plates and shells containing through cracks may be found in 
[21], [22], [24]-[26], For example, if we use a length parameter a* (usually 
a half crack length) to normalize the dimensions, coordinates, and other 
quantities as 

x-j/a* = x, u 2 /a* = v, a p /a* = a p , b p /a* = b p (44) 

and define 

8 2 (x,+0) = ^(x) , v(x,+0) = <f> 2 (x) , (45) 

The integral equations (35) for an infinite plate may be expressed as [31] 



P , P 2 hr} 1 

z j ^ 2 j( x i »t)fj p (t)dt - -g- c 2-j (x-j ) j f i r (t )dt 


c 22 ( X l) 


f 2 r(t)dt ~ £ » (a r <x.| <b r , r 1,...,P) , 


rr 


f ir (t)dt = 0, (1=1,2, r=l,2,...,P) 
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4>p ( 5 ) 2 P , ■ 

— - — de = ^ [-o«+a(x)L x e E (a p ,b p ) , 


( 47 ) 


I 



where y = 10(a*) 2 /h 2 . 

To give an idea about the applications of the line spring model some 
sample results are shown in Tables 5-10 and in Figures 6-13. Table 5 gives 
the normalized stress intensity factor at the deepest penetration points of 
two coplanar semielliptic surface cracks symmetrically located on the oppo- 
site sides of an infinite plate under uniform tension a Q =N/h (Fig. 4). The 
table also shows the corresponding plane strain result which is the limiting 
value of the stress intensity factor as the crack length 2a tends to infinity. 
In this as well as in other examples discussed in this section the crack 
border in x-^ plane is assumed to be defined by (Fig. 1) 


L( X] ) = L 0 /T^xf/iF (48) 

where a is the half crack length. 

The results for the two semi-elliptic coplanar surface cracks of same 
dimensions and located on the same side of an infinite plate are given in 
Tables 6 and 7. Table 6 shows the maximum value of the normalized stress 
intensity factor for the plate under uniform tension o Q in a direction per- 
pendicular to the plane of the crack. The notation used in these tables 
regarding the dimensions is the same as in Fig. 5a with d=h/2. The stress 
intensity factor for the same plate under bending is given in Table 7. The 
tables also show the stress intensity factors for the single surface crack 
which are the limiting values of the two crack results for b-«°. In the 
two crack problem considered the distribution of stress intensity factor 
K(x-j) along the crack border (or as a function of x-j ) is somewhat skewed 
and the maximum K occurs at a value of x-| that is somewhat less than b+a 
[24]. 

Figure 6 shows the normalized stress intensity factor in a plate con- 
taining a symmetrically located internal elliptic crack (d=0. Fig. 3) at 
the points where the minor axis intersects the crack border. The results 
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are given for a plate under uniform tension perpendicular to the plane of 
the crack. In this as well as in the symmetric surface crack problem con- 
sidered in Table 5 there is no bending and hence the solution by the line 
spring method is extremely simple. The figure also shows the finite element 
results from [3], Disregarding some small values of a/L Q for which the 
line spring model is not really suitable, it may be noted that the agreement 
between the two results is fairly good. 

Some sample results for an excentrically located internal elliptic crack 
are shown in Figures 7-10. K^ and Kg shown in these figures correspond to 
the points at the midsection of the ellipse (see the insert in Fig. 7). 

Figures 7 and 8 show the normalized stress intensity factors in a plate 
under uniform tension a Q for a fixed value of a/L Q =4 and varying values of 
d/h and L Q /h. Figures 9 and 10 show the tension results for d/h = 0.15 and 
varying values of L Q /h and a/L Q . 

The normalized stress intensity factor for two symmetrically located 
identical coplanar internal elliptic cracks in a plate under uniform tension 
is given in Table 8 (d=0. Fig. 5a). The table shows the results at the 
midsection of the ellipses. The result for the limiting case of a single 
crack is also given in the table. 

The results for three identical cracks shown in Fig. 5b are given in 
Table 9. In this case, too, K^ and Kg refer to the points of intersection 
of the minor axes of ellipses with the crack border (Fig. 5b), Since d^O 
K^ and Kg are not equal. One may observe that the stress intensity factors 
for the middle crack are only slightly higher than that for the two end 
cracks. If one may make one general observation regarding the interaction 
between multiple cracks, it would be that for the same crack lengths and 
distances in x-j direction, the interaction for the part-through cracks is 
much weaker than the interaction between through cracks. 

An example for the distribution of the stress intensity factors along 
the crack border is given in Table 10. More extensive results on the inter- 
action of multiple part-through cracks of various geometries in an infinite 
plate may be found in [30], 

Some results for a plate of finite width are shown in Figures 11 and 
12. Fig. 11 shows the distribution of the normalized stress intensity factors 
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for a symmetrically located surface crack having a semi-elliptic or a rec- 
tangular profile. The normalizing stress intensity factors shown in these 
figures are the corresponding plane strain values for an edge-cracked strip 
and are defined by (see equations 8 and 9) 


to = Hr ^ 9 t (s o> > K I 


6M 


bo 


'W' ^ W’ s o = L o /h 


(49) 


Figure 12 shows an example for the corner cracks. More extensive results 
for multiple cracks in a plate of finite width obtained by using the line 
.spring model and comparison with some of the finite element solutions may 
be found in [24]. 

Aside from some additional rather complicated Fredholm kernels in the 
integral equations for shells, from a viewpoint of applications of the line 
spring model, the problems in shells and plates are identical. The results 
for various crack and shell geometries obtained by using a transverse shear 
theory of shallow shells are given in [22], [23], [25] and [26]. Figure 13 
shows an example for a pressurized cylindrical shell with a fixed end con- 
taining a semi-elliptic axial surface crack. The problem may simulate a 
rigid end plate or a relatively heavy flange. It is seen that the effect 
of shell curvature on the stress intensity factors can be very significant. 

In reviewing the results one may conclude that, despite its simplicity, 
carefully and judiciously applied, the line spring model may give very use- 
ful results for some three-dimensional part-through crack problems that are 
otherwise analytically intractable. The method is naturally suited to 
account for plastic deformations in certain approximate ways. The questions 
currently being studied concern the extension of the model to mixed-mode 
problems where, unlike the Mode I case. Modes II and III are always coupled. 
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Table 1 


. The stress intensity factors and the coefficients bj and Cj of 
the shape functions g t (s) and g^s) in a strip containing an 
edge crack of length L and subjected to uniform tensile stress 
N/h and bending moment M (see eqs. 8 and 9 and Fig. 2). 


Tension Bending 


L/h 

K-|/(N/h)/nT 

K-j / ( 6M/ h 2 ) /ttU 

b i 

c i 

i 

-> 0 

1.12152226 

1 .12152226 

1.1215 

1 .1202 

0 

10" 5 

1.121522 

1.1215 


-1.8872 

1 

10- 3 

1.121531 

1 .1202 

-12.3877 

18.0143 

2 

0.1 


1.0472 


-87.3851 

3 

0.2 

1.3673 


-188.6080 

241.9124 

4 

0.3 

1.6599 

1.1241 

207.3870 

-319.9402 

5 

0.4 

2.1114 

1.2606 


168.0105 1 

6 

0.5 

2.8246 

1.4972 




0.6 


1.9140 




0.7 

6.3549 





0.8 

11.955 

4.6764 




0.85 


6.9817 




0.9 


12.462 




0.95 

99.14 

34.31 





Table 2. The stress intensity factors and the coefficients of shape 
functions in a strip with a symmetric internal crack (d=0. 
Fig. 3b) under tension and bending (see eqs. 27 and 28). 


Tension Bending 


L/h 

K^/ o/7tL/2 

K A /m MTZ 

b Ai =b Bi 

c Ai c Bi 

i 

0.05 

1.0002 

0.0250 

0.7071 

0.1013 

0 

0.1 

1.0060 

0.0500 

0.4325 

-0.4629 

1 

0.2 

1.0246 

0.1001 

-0.1091 

15.0622 

2 

0.3 

1.0577 

0.1505 

7.3711 

-143.7384 

3 

0.4 

1.1094 

0.2023 

-57.7894 

807.2449 

4 

0.5 

1.1867 

0.2573 

271.1551 

-2844.8525 

5 

0.6 

1.3033 

0.3197 

-744.4204 

6468.9152 

6 

0.7 

1.4884 

0.3986 

1183.9529 

-9477.5512 

7 

0.8 

1.8169 

0.51 86 

-1001.4920 

8638.7826 

8 

0.9 

2.585 

0.7776 

347.9786 

-4455.2167 

9 

0.95 

4.252 

1.1421 

347.9786 

5959.4888 

10 


20 - 







Table 3. Coefficients of the shape functions in a strip with an asymmetric internal 
crack under tension and bending; d^O, see eqs. (27) and (28) and Fib. 3b. 


d/h 

i 

b Ai 

b Bi 

c Ai 

c Bi 


0 

0.7071 

0.7071 

0.0708 

0.0707 


1 

0.4597 

0.4347 

-0.0623 

-0.3701 


2 

0.7671 

-0.0915 

13.1229 

0.5654 


3 

0.1552 

2.6973 

-166.4280 

-6.6423 


4 

-9.3017 

-14.1195 

1145.8217 

45.7189 

0.05 

5 

97.3172 

54.9653 

-4762.0914 

-189.9515 

6 

-413.9673 

-135.3432 

12511.5152 

498.8463 


7 

936.4719 

205.3051 

-20927.0019 

-834.5704 


8 

-1078.2322 

-173.3480 

21613.9362 

862.1672 


9 

504.0555 

62.8847 

-12568.0268 

-501.4354 





3148.4879 

125.5869 


0 

0.7071 

0.7072 

0.1415 

0.1414 


1 

0.5498 

0.5043 

-0.1734 

-0.3871 



1.5235 

-0.5779 

18.7434 

1.2936 



-2.2395 

7.6480 

-266.7713 

-17.0715 



-5.2844 

-52.8793 

2066.4692 

132.0282 

0.1 


226.0267 

257.2074 

-9661.5218 

-617.3023 



-1423.2887 

-799.7410 

28556.2764 

1826.3191 



4348.1446 

1530.8314 

-53734.1216 

-3441.9797 



-6553.5540 

-1634.0240 

62435.9340 

4007.6642 



3959.2116 

749.0673 

-40844.2364 ; 

-2628.6642 


10 



11511,5912 

743.3343 


0 

0.7071 

0.7072 

0.2122 

0.2121 


1 

0.7028 

0.6376 

-0.2929 

-0.4042 


2 

2.7653 

-1.2331 

26.3239 

2.2494 


3 

-7.2036 

19.0057 

-427.2558 

-33.6757 


4 

9.1384 

-173.8407 

3782.9591 

297.4990 

0.15 

5 

667.4954 

1108.9410 

-20214.1250 

-1590.1109 

6 

-6105.7233 

-4517.1019 

68285.5344 

5378.6049 


7 

25260.2847 

11317.3469 

-146859.5866 

-11588.2217 


8 

-50586.0954 

-15802.5485 

195038.2341 

15425.1699 


9 

40325.8388 

9475,7480 

-145833.6228 

-11566.8288 


10 



46980.5243 

3740.3538 
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Table 3 - cont. 


d/h 

i 

b Ai 

b Bi 

C Ai 

c Bi 


0 

0.7071 

0.7072 

0.2829 

0.2828 


1 

| 0.9394 

0.8534 

-0.4105 

-0.4192 


2 

5.0186 

-2.2518 

36.3675 

3.4524 


3 

-19.6345 

47.2610 

-686.9924 

-59.3848 


4 

76.1489 

-589.3736 

7097.1745 

611.7194 

0.20 

5 

2376.8770 

5125.5432 

-44245.1037 

-3814.7743 


6 

-32402.0663 

-28413.7181 

174386.6733 

15058.0019 


7 

187563.3073 

96818.0664 

-437594.1661 

-37860.2290 


8 

-517758.7465 

-183743.9141 

678087.6506 

58816.9514 


9 

566112.6482 

149736.5141 

-591607.7634 

-51479.6217 


10 



222394.8277 

19433.7456 


0 

0.7071 

0.7071 

0.3536 

0.3536 


1 

1 .3041 

1.1918 

-0.5246 

-0.4284 


2 

9.6510 

-4.1793 

50,2589 

4,8592 


3 

-57.8163 

129.7358 

-1136.1027 

-97.6558 


4 

450.7761 

-2325.3551 

14085.0245 

1206.5210 

0.25 

5 

10351.8621 

29069.6053 

-105370.0173 

-9029.9162 


6 

-227973.9193 

-231101.0011 

498386,1475 

42788.3389 


7 

1960136.6080 

1128028.3541 

-1500844.1305 

-129153.3052 


8 

-7899243.5583 

-3063223.0225 

2791092.0497 

240895,9775 


9 

12538854.7550 

3568834.5097 

-2922538,2779 

-253169,5873 


10 



1318591.9416 ! 

114779.2965 


0 

0.7071 

0.7071 

0.4244 

0.4243 


1 

1.9027 

1.7480 

-0.6553 

-0.4321 


2 

20.8636 

-8.9087 

72.3179 

6.6403 


3 

-197.8895 

440.7103 

-2037.1553 

-160.0429 


4 

2675.1513 

-12390.9668 

31569.7505 

2468.7747 

0.30 

5 

71601.3880 

241718.2677 

-295202.5115 

-23090.8486 

6 

-2639273.5314 

-2989605.6659 

1745346.7908 

136788.4159 


7 

35994016.1511 

22618755.1306 

-6570139.7052 

-516213.4074 


8 

-227958779.2155 

-94950523.9757 

15273818.2302 

1203928.2026 


9 

567162515.8604 

170542044.2578 

-19993047.6604 

-1582260.0614 


10 



11276922.9053 

897240.7983 
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Table 4, The stress intensity factor and the coefficients 

of the shape function for a strip with two coll inear 
symmetric edge cracks of depths L/2 under tension a. 


L/h 


bi 

i 

0.0001 

1.1221 

0.7934 

0 

0.1 

1,1231 

0.0775 

1 

0.2 

1.1254 

-0.7542 

2 

0.3 

1.1292 

7.5825 

3 

0.4 

1.1370 

-12.1712 

4 

0.5 

1.1546 

-186.5011 

5 

0.6 

1.2117 

1236.2858 

6 

0.7 

1.3254 

-3043.6190 

7 

0.8 

1.5393 

3350.3456 

8 

0.9 

2.0836 

-1374.8426 

9 


Table 5. Normalized stress intensity factor K/K 0 calculated at the midsection 



of two opposite planar 

" ■ f ■ 

elliptic surface 

cracks in 

a plate under uni 


form tension 

ff o» K o =cr o 

/iFy? 

(Fig. 4). 




a/L Q -*■ 1.0 

2.0 

3.0 

4.0 

10.0 

100.0 

Plane 

V h 







Strain 

0.1 

1.060 

1.089 

1.099 

1.104 

1.113 

1.119 

1.123 

0.2 

1.009 

1.062 

1.081 

1.091 

1.109 

1.121 

1.125 

0.3 

0.966 

1.028 

1.065 

1.079 

1.106 

1.124 

1.129 

0.4 

0.929 

1.019 

1.053 

1.072 

1.108 

1.131 

1 .137 

0.5 

0.902 

1.008 

1.050 

1.073 

1.118 

1.148 

1.155 

0.6 

0.902 

1.038 

1.080 

1.108 

1.165 

1.203 

1.212 

0.7 

0.929 

1.082 

1.149 

1.186 

1.262 

1.315 

1.325 

0.8 

0.997 

1.195 

1.284 

1.336 

1.445 

1.524 

1.539 
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Table 6. Maximum normalized stress intensity factor K/K 0 for two 
planar elliptic su rfac e cracks in a plate under uniform 
tension a Q , K^a /tt’L . (Fig. 5a, d=h/2) d=h/2 



0.1 2 0.996 
0.1 4 1.072 
0.1 10 1.127 


0.2 2 0.979 
0.2 4 1.107 
0.2 10 1.221 

0.3 2 1.007 
0.3 4 1.189 
0.3 10 1.382 

0.4 2 1.051 
0.4 4 1.299 
0.4 10 1.600 

0.5 2 1.106 
0.5 4 1.430 
0.5 10 1.879 

0.6 2 1.230 
0.6 4 1.568 
0.6 10 2.208 

0.7 2 1.370 
0.7 4 1.694 
0.7 10 2.532 


1.0 

4.0 

20.0 

0.986 

0.982 

0.981 

1.065 

1.062 

1.062 

1.122 

1.121 

1.120 

0.959 

0.951 

0.949 

1.090 

1.084 

1.082 

1.210 

1.207 

1.206 

0.975 

0.963 

0.961 

1.160 

1.151 

1.149 

1.364 

1.359 

1.358 

1.006 

0.991 

0.989 

1.255 

1.243 

1.240 

1.571 

1.563 

1.562 

1.048 

1 .032 

1.030 

1.370 

1.354 

1.352 

1 .836 

1.824 

1.822 

1.096 

1 .080 

1.077 

1,491 

1.472 

1.469 

2.141 

2,124 

2.121 

1.120 

1.100 

1,097 

1.575 

1.553 

1.550 

2.434 

2.409 

2.405 


Single 

Crack 


0.981 

1.062 

1.120 

0.949 

1.082 

1.206 

0,961 

1.149 

1.358 

0.989 

1.240 

1,562 

1.030 

1.352 

1.821 

1.077 

1.469 

2.121 

1,097 

1.550 

2.405 
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Table 7 


. Maximum normalized stress intensity factor K/Kq for two 
planar elliptic s urfa ce cracks in a plate under uniform 
bending M, I<o=a b AT^, a b =6M/h 2 (Fig, 5a, d=h/2). 


L o /h 

by 

a/L o 

'a -v 0.1 

1.0 

4.0 

20.0 

Single 

Crack 

0.1 

2 

0.874 

0.864 

0.861 

0.860 

0,860 

0.1 

4 

0.943 

0.936 

0.934 

0.933 

0.933 

0.1 

10 

0.992 

0.988 

0.987 

0.986 

0.986 

0.2 

2 

0.766 

0.728 

0.719 

0.718 

0.718 

0.2 

4 

0.847 

0.830 

0.825 

0.824 

0.824 

0.2 

10 

0.940 

0.930 

0.927 

0.927 

0.927 

0.3 

2 

0.751 

0.665 

0.651 

0.650 

0.650 

0.3 

4 

0.803 

0.755 

0.745 

0.744 

0.744 

0.3 

10 

0.915 

0.900 

0.896 

0.896 

0.895 

0.4 

2 

0.792 

0.677 

0.659 

0.658 

0.656 

0.4 

4 

0.801 

0.726 

0.713 

0.711 

0.711 

0.4 

10 

0.923 

0.902 

0.896 

0.895 

0.895 

0.5 

2 

0.826 

0.684 

0,663 

0.661 

0.659 

0.5 

4 

0.834 

0.719 

0.703 

0.701 

0.700 

0.5 

10 

0.950 

0.910 

0.902 

0.901 

0.901 

0.6 

2 

0.855 

0.686 

0.662 

0.660 

0,659 

0.6 

4 

0.909 

0.743 

0.724 

0.722 

0.721 

0.6 

10 

0.995 

0.925 

0.912 

0.910 

0.910 

0.7 

2 

0.874 

0.683 

0.658 

0.655 

0,654 

0.7 

4 

0.989 

0.784 

0.761 

0.759 

0.757 

0.7 

10 

1.064 

0.956 

0.939 

0,937 

0.936 
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Table 8. 


Normalized stress intensity factor K/K 0 calculated at the 
midsection of symmetrically located (d=0) two identical 



planar internal 

elliptic 

cracks in 

a plate 

under uniform 


tension a_, K = 
O’ o 

a./mL /2 
0 0 

(Fig. 5a). 



b/a -+ 0.1 

1.0 

4.0 

20.0 

Single 

L o /h 

a ' L 0 





Crack 

0.1 

2 

0.977 

0.976 

0.976 

0.976 

0.976 

0.1 

4 

0.987 

0.987 

0.987 

0.986 

0.987 

0.1 

10 

0.994 

0.993 

0.993 

0.993 

0.993 

0.2 

2 

0.975 

0.972 

0.971 

0.971 

0.971 

0.2 

4 

0.995 

0.994 

0.993 

0.993 

0.993 

0.2 

10 

1.008 

1.007 

1.007 

1 .007 

1.007 

0.3 

2 

0.980 

0.980 

0.979 

0.979 

0.979 

0.3 

4 

1.015 

1 .013 

1.012 

1.012 

1.012 

0.3 

10 

1.035 

1.034 

1 .034 

1.034 

1.034 

0.4 

2 

1.007 

1.002 

1.000 

1.000 

1.000 

0.4 

4 

1.051 

1.048 

1.047 

1.047 

1.047 

0.4 

10 

1.080 

1.078 

1.078 

1.078 

1.078 

0.5 

2 

1.046 

1.039 

1.037 

1.036 

1.036 

0.5 

4 

1.106 

1.102 

1.101 

1.101 

1.101 

0.5 

10 

1.147 

1.145 

1.145 

1.145 

1.145 

0.6 

2 

1.108 

1.098 

1.095 

1.095 

1,095 

0.6 

4 

1.190 

1.185 

1.183 

1.183 

1.183 

0.6 

10 

1.248 

1 .246 

1.245 

1.245 

1.245 

0.7 

2 

1 1.205 

1.192 

1.188 

1.187 

1.187 

0.7 

4 

1.321 

1.313 

1.311 

1.310 

1.310 

0.7 

10 

1.407 

1.404 

1.403 

1.403 

1.403 

0.8 

2 

1.367 

1.348 

1.342 

1.341 

1.341 

0.8 

4 

1.541 

1.529 

1.526 

1.525 

1.525 

0.8 

10 

1.681 

1.676 

1,674 

1.674 

1.674 

0.9 

2 

1.703 

1.672 

1.662 

1.661 

1.661 

0.9 

4 

2.007 

1.988 

1.982 

1.981 

1.981 

0.9 

io ; 

2.285 

2.275 

2.272 

2.272 



-26 


Table 9. Normalized stress intensity factors K^/Kq and Kg/K 0 calculated 
at the midsection of three identical planar internal elliptic 
cracks in a plate under uniform tension a 0 , Kq^q^VL. /2 . 

The middle crack 


b/a -> 0.1 1.0 4.0 20.0 



d 

L o /h 

a/L o 

• 1 

1 • U 

"t • u 

C.KJ • U 


0.1 

0.1 

2 

0.983 

0.981 

0.979 

0.979 


0.1 

0.1 

4 

0.994 

0.992 

0.992 

0.991 


0.1 


2 

0.985 

0.980 

0.977 

0.976 


0.1 

ESS 

4 

1.007 

1.004 

1.002 

1.001 


0.1 

0.3 

2 

, 1.006 

0.997 

0.991 

0.989 

k a 

0.1 

0.3 

4 

1.040 

1.035 

1,032 

1.031 

K o 

0.2 

0.1 

2 

0.981 

0.978 

0.976 

0.975 


0.2 

0.1 

4 

0.995 

0.994 

0,992 

0,992 


0.2 

0.2 

2 

0.993 

0.986 

0.981 

0,979 


0.2 

0.2 

4 

1,023 

1.018 

1.015 

1.014 


0.3 

0.1 

2 i 

0.983 

0.979 

0.976 

0.975 


0.3 

0.1 

4 

1.003 

1.001 

0.999 

0,998 


0.1 

0.1 

2 

0.983 

0.981 

0.979 

0.979 


0.1 

0.1 

4 

0.994 

0.992 

0.991 

0.991 


0.1 

0.2 

i 

2 

0.984 

0.980 

0.977 

0.976 


0.1 

0.2 

4 

1,003 

1.001 

0.999 

0,999 


0.1 

0.3 

2 

0.997 

0.990 

0.986 

0.984 


0.1 

0.3 

4 

1.026 

1.022 

1.020 

1,019 


0.2 

0.1 

2 

0.981 

0.978 

0.976 

0,975 

k b 

0.2 

0.1 

4 

0.994 

0.993 

0.992 

0,991 


0.2 

0.2 

2 

0.987 

0.981 

0.977 

0.976 


0.2 

0.2 

4 

1.012 

1.009 

1.006 

1.005 


0.3 

0.1 

2 

0.981 

0.977 

0.975 

0,974 


0.3 

0.1 

4 

0.999 

0.997 

0.996 

0.995 
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Table 9 - Cont 


The outer cracks 


b/a -v 0.1 1.0 4.0 20.0 



u 

L o /M 

d / L 0 



0.1 

0.1 

' 2 


0.980 

0.979 

0.979 


0.1 

0.1 

j 4 

1 0.993 


0.991 

0.991 


0.1 

0.2 

2 

0.982 

0.978 

0.976 

0.976 


0.1 

0.2 

4 

1.005 

1.003 

1.002 

1.001 


0.1 

i 0.3 

2 

0.999 

0.994 

0.991 

0.989 

k a 

0.1 

0.3 

4 

1.037 

1.034 

1.032 

1.031 

K 

0 

0.2 

0.1 

2 

0.979 

0.977 

0.976 

0.975 


0.2 

0.1 

4 

0.994 

0.993 

0.992 

0.992 


0.2 

0.2 

2 

0.988 

0.983 

0.980 

0.979 


0.2 

0.2 

4 

1.019 

1.016 

1.014 

1.014 


0.3 

0.1 

2 

0.980 

0.978 

0.976 

0.975 


0.3 

0.1 

4 

1.001 

1.000 

0.999 

0.998 


0.1 

0.1 

2 

0.982 

0.980 

0.979 

0.979 


0.1 

0.1 

4 

0.993 

0.992 

0.991 

0.991 


0.1 

0.2 

2 

0.981 

0.978 

0.976 

0.976 


0.1 

0.2 

4 

1.002 

1.000 

0.999 

0.999 


0.1 

0.3 

2 

0.992 

0.988 

0.985 

0.984 

k b 

0.1 

0.3 

4 

1.023 

1.021 

1.020 

1.019 

K n ' 

0.2 

0.1 

2 

0.979 

0.978 

0.976 

0.975 

0 

0.2 

0.1 

4 

0.993 

0.992 

0.991 

0.991 


0.2 

0.2 

2 

0.983 

0.979 

0.977 

0.976 


0.2 

0.2 

4 

1.010 

1.007 

1.006 

1.005 


0.3 

0.1 

2 

0.978 

0.976 

0.974 

0.974 


0.3 

0.1 

4 

. 

0.998 

0.996 

0.995 

0.995 
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Table 10 


, Normalized stress intensity factors on the crack 
front for an internal elliptic crack in a plate 
under uniform membrane load Np2=ho 0 and bending 
moment M?p=h 2 o h/h with d/h=0 .20,. L o /h=0.45, a/L 0 = 
4, K=a 0 Al /2 or K =a b /ift /2 (see Fig. 3), 


x-|/a 

<Wn 

Wn 

( Wm 

( Wm 

0.90 

0.723 

0.707 

0.357 

0.218 

0.80 

0.865 

0.831 

0.454 

0.229 

0.70 

0.974 

0.916 

0.531 

0.233 

0.60 

1.065 

0.980 

0.595 

0.233 

0.50 

1.143 

1.028 

0.649 

0.233 

0.40 

1.211 

1.068 

0.697 

0.233 

0.30 

1.266 

1.096 

0.735 

0,233 

0.20 

1.307 

1.116 

0.763 

0,233 

0.10 

1.333 

1.128 

0.781 

0,233 

0.00 

1.342 

1.132 

0.787 

0.232 
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e 



2 The corresponding plane strain problem 
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Fig. 4 Geometry and notation for two coplanar semi-elliptic surface 
cracks on the opposite sides of the plate 
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Fig. 5 Geometry and notation for.coplanar internal elliptic cracks 
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L 0 /h 


Fig. 7 Normalized stress intensity factor at the midsection of an 

eccentrically located internal ell iptic crack in a plate under 
uniform tension (a/L 0 =4, K 0 =a 0 v / TrL 0 /2 , subscript A refers 
to the smaller net ligament). 
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U/h =0.6 



Stress intensity factor at the midsection of a semi-elliptic 
axial surface crack in a pressurized cylinder with a fixed 
end (v=0.3, a/h=l , L o /h=0,4) 
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